Phenotypic trait variation in a long-term multisite common garden experiment of Scots pine in Scotland.

Multisite common garden experiments, exposing common pools of genetic diversity to a range of environments, allow quantification of plastic and genetic components of trait variation. For tree species, such studies must be long term as they typically only express mature traits after many years. As well as evaluating standing genetic diversity, these experiments provide an ongoing test of genetic variation against changing environmental conditions and form a vital resource for understanding how species respond to abiotic and biotic variation. Finally, quantitative assessments of phenotypic variation are essential to pair with rapidly accumulating genomic data to advance understanding of the genetic basis of trait variation, and its interaction with climatic change. We describe a multisite, population-progeny, common garden experiment of the economically and ecologically important tree species, Scots pine, collected from across its native range in Scotland and grown in three contrasting environments. Phenotypic traits, including height, stem diameter and budburst were measured over 14 growing seasons from nursery to field site. The datasets presented have a wide range of applications.

Identifying and testing marker-trait associations for growth and phenology in three pine species: Implications for genomic prediction.

In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 single nucleotide polymorphisms (SNPs) from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments.

Taming the massive genome of Scots pine with PiSy50k, a new genotyping array for conifer research.

Pinus sylvestris (Scots pine) is the most widespread coniferous tree in the boreal forests of Eurasia, with major economic and ecological importance. However, its large and repetitive genome presents a challenge for conducting genome-wide analyses such as association studies, genetic mapping and genomic selection. We present a new 50K single-nucleotide polymorphism (SNP) genotyping array for Scots pine research, breeding and other applications. To select the SNP set, we first genotyped 480 Scots pine samples on a 407 540 SNP screening array and identified 47 712 high-quality SNPs for the final array (called 'PiSy50k'). Here, we provide details of the design and testing, as well as allele frequency estimates from the discovery panel, functional annotation, tissue-specific expression patterns and expression level information for the SNPs or corresponding genes, when available. We validated the performance of the PiSy50k array using samples from Finland and Scotland. Overall, 39 678 (83.2%) SNPs showed low error rates (mean = 0.9%). Relatedness estimates based on array genotypes were consistent with the expected pedigrees, and the level of Mendelian error was negligible. In addition, array genotypes successfully discriminate between Scots pine populations of Finnish and Scottish origins. The PiSy50k SNP array will be a valuable tool for a wide variety of future genetic studies and forestry applications.

Candidate Genes for the High-Altitude Adaptations of Two Mountain Pine Taxa.

Mountain plants, challenged by vegetation time contractions and dynamic changes in environmental conditions, developed adaptations that help them to balance their growth, reproduction, survival, and regeneration. However, knowledge regarding the genetic basis of species adaptation to higher altitudes remain scarce for most plant species. Here, we attempted to identify such corresponding genomic regions of high evolutionary importance in two closely related European pines, Pinus mugo and P. uncinata, contrasting them with a reference lowland relative-P. sylvestris. We genotyped 438 samples at thousands of single nucleotide polymorphism (SNP) markers, tested their genetic differentiation and population structure followed by outlier detection and gene ontology annotations. Markers clearly differentiated the species and uncovered patterns of population structure in two of them. In P. uncinata three Pyrenean sites were grouped together, while two outlying populations constituted a separate cluster. In P. sylvestris, Spanish population appeared distinct from the remaining four European sites. Between mountain pines and the reference species, 35 candidate genes for altitude-dependent selection were identified, including such encoding proteins responsible for photosynthesis, photorespiration and cell redox homeostasis, regulation of transcription, and mRNA processing. In comparison between two mountain pines, 75 outlier SNPs were found in proteins involved mainly in the gene expression and metabolism.

Development of a single nucleotide polymorphism array for population genomic studies in four European pine species.

Pines are some of the most ecologically and economically important tree species in the world, and many have enormous natural distributions or have been extensively planted. However, a lack of rapid genotyping capability is hampering progress in understanding the molecular basis of genetic variation in these species. Here, we deliver an efficient tool for genotyping thousands of single nucleotide polymorphism (SNP) markers across the genome that can be applied to genetic studies in pines. Polymorphisms from resequenced candidate genes and transcriptome sequences of P. sylvestris, P. mugo, P. uncinata, P. uliginosa and P. radiata were used to design a 49,829 SNP array (Axiom_PineGAP, Thermo Fisher). Over a third (34.68%) of the unigenes identified from the P. sylvestris transcriptome were represented on the array, which was used to screen samples of four pine species. The conversion rate for the array on all samples was 42% (N = 20,795 SNPs) and was similar for SNPs sourced from resequenced candidate gene and transcriptome sequences. The broad representation of gene ontology terms by unigenes containing converted SNPs reflected their coverage across the full transcriptome. Over a quarter of successfully converted SNPs were polymorphic among all species, and the data were successful in discriminating among the species and some individual populations. The SNP array provides a valuable new tool to advance genetic studies in these species and demonstrates the effectiveness of the technology for rapid genotyping in species with large and complex genomes.

Invasion, isolation and evolution shape population genetic structure in Campanula rotundifolia.

The distribution and genetic structure of most plant species in Britain and Ireland bear the imprint of the last ice age. These patterns were largely shaped by random processes during recolonization but, in angiosperms, whole-genome duplication may also have been important. We investigate the distribution of cytotypes of Campanula rotundifolia, considering DNA variation, postglacial colonization, environmental partitioning and reproductive barriers. Cytotypes and genome size variation from across the species' range were determined by flow cytometry and genetic variation was assessed using cpDNA markers. A common garden study examined growth and flowering phenology of tetraploid, pentaploid and hexaploid cytotypes and simulated a contact zone for investigation of reproductive barriers. Irish populations were entirely hexaploid. In Britain, hexaploids occurred mostly in western coastal populations which were allopatric with tetraploids, and in occasional sympatric inland populations. Chloroplast markers resolved distinct genetic groups, related to cytotype and geographically segregated; allopatric hexaploids were distinct from tetraploids, whereas sympatric hexaploids were not. Genome downsizing occurred between cytotypes. Progeny of open-pollinated clones from the contact zone showed that maternal tetraploids rarely produced progeny of other cytotypes, whereas the progeny of maternal hexaploids varied, with frequent pentaploids and aneuploids. The presence of distinctive hexaploid chloroplast types in Ireland, Scottish islands and western mainland Britain indicates that its establishment preceded separation of these land masses by sea-level rise c. 16 000 years BP. This group did not originate from British tetraploids and probably diverged before postglacial invasion from mainland Europe. The combination of cytotype, molecular, contact zone and common garden data shows an overall pattern reflecting postglacial colonization events, now maintained by geographic separation, together with more recent occasional local in situ polyploidisation. Reproductive barriers favour the persistence of the tetraploid to the detriment of the hexaploid.

Early phenology and growth trait variation in closely related European pine species.

Closely related taxa occupying different environments are valuable systems for studying evolution. In this study, we examined differences in early phenology (bud set, bud burst) and early growth in a common garden trial of closely related pine species: Pinus sylvestris, P. mugo, and P. uncinata. Seeds for the trial were sourced from populations across the ranges of each species in Europe. Over first 4 years of development, clear differences were observed between species, while the most significant intraspecific differentiation was observed among plants from P. sylvestris populations from continental European locations. Trait differences within P. sylvestris were highly correlated with altitude and latitude of the site of origin. Meanwhile, P. mugo populations from the Carpathians had the earliest bud set and bud flush compared to other populations of the species. Overall, populations from the P. mugo complex from heterogeneous mountain environments and P. sylvestris from the Scottish Highlands showed the highest within-population variation for the focal traits. Although the three species have been shown to be genetically highly similar, this study reveals large differences in key adaptive traits both among and within species.

Molecular signatures of divergence and selection in closely related pine taxa.

Efforts to detect loci under selection in plants have mostly focussed on single species. However, assuming that intraspecific divergence may lead to speciation, comparisons of genetic variation within and among recently diverged taxa can help to locate such genes. In this study, coalescent and outlier detection methods were used to assess nucleotide polymorphism and divergence at 79 nuclear gene fragments (1212 SNPs) in 16 populations (153 individuals) of the closely related, but phenotypically and ecologically distinct, pine taxa Pinus mugo, P. uliginosa and P. uncinata across their European distributions. Simultaneously, mitochondrial DNA markers, which are maternally inherited in pines and distributed by seeds at short geographic distance, were used to assess genetic relationships of the focal populations and taxa. The majority of nuclear loci showed homogenous patterns of variation between the taxa due to a high number of shared SNPs and haplotypes, similar levels of polymorphism, and low net divergence. However, against this common genetic background and an overall low population structure within taxa at mitochondrial markers, we identified several genes showing signatures of selection, accompanied by significant intra- and interspecific divergence. Our results indicate that loci involved in species divergence may be involved in intraspecific local adaptation.

Reconstructing the plant mitochondrial genome for marker discovery: a case study using Pinus.

Whole-genome-shotgun (WGS) sequencing of total genomic DNA was used to recover ~1 Mbp of novel mitochondrial (mtDNA) sequence from Pinus sylvestris (L.) and three members of the closely related Pinus mugo species complex. DNA was extracted from megagametophyte tissue from six mother trees from locations across Europe, and 100-bp paired-end sequencing was performed on the Illumina HiSeq platform. Candidate mtDNA sequences were identified by their size and coverage characteristics, and by comparison with published plant mitochondrial genomes. Novel variants were identified, and primers targeting these loci were trialled on a set of 28 individuals from across Europe. In total, 31 SNP loci were successfully resequenced, characterizing 15 unique haplotypes. This approach offers a cost-effective means of developing marker resources for mitochondrial genomes in other plant species where reference sequences are unavailable.

Genetic structure in the European endemic seabird, Phalacrocorax aristotelis, shaped by a complex interaction of historical and contemporary, physical and nonphysical drivers.

Geographically separated populations tend to be less connected by gene flow, as a result of physical or nonphysical barriers preventing dispersal, and this can lead to genetic structure. In this context, highly mobile organisms such as seabirds are interesting because the small effect of physical barriers means nonphysical ones may be relatively more important. Here, we use microsatellite and mitochondrial data to explore the genetic structure and phylogeography of Atlantic and Mediterranean populations of a European endemic seabird, the European shag, Phalacrocorax aristotelis, and identify the primary drivers of their diversification. Analyses of mitochondrial markers revealed three phylogenetic lineages grouping the North Atlantic, Spanish/Corsican and eastern Mediterranean populations, apparently arising from fragmentation during the Pleistocene followed by range expansion. These traces of historical fragmentation were also evident in the genetic structure estimated by microsatellite markers, despite significant contemporary gene flow among adjacent populations. Stronger genetic structure, probably promoted by landscape, philopatry and local adaptation, was found among distant populations and those separated by physical and ecological barriers. This study highlights the enduring effect of Pleistocene climatic changes on shag populations, especially within the Mediterranean Basin, and suggests a role for cryptic northern refugia, as well as known southern refugia, on the genetic structure of European seabirds. Finally, it outlines how contemporary ecological barriers and behavioural traits may maintain population divergence, despite long-distance dispersal triggered by extreme environmental conditions (e.g. population crashes).

Has Scots pine (Pinus sylvestris) co-evolved with Dothistroma septosporum in Scotland? Evidence for spatial heterogeneity in the susceptibility of native provenances.

Spatial heterogeneity in pathogen pressure leads to genetic variation in, and evolution of, disease-related traits among host populations. In contrast, hosts are expected to be highly susceptible to exotic pathogens as there has been no evolution of defence responses. Host response to pathogens can therefore be an indicator of a novel or endemic pathosystem. Currently, the most significant threat to native British Scots pine (Pinus sylvestris) forests is Dothistroma needle blight (DNB) caused by the foliar pathogen Dothistroma septosporum which is presumed to be exotic. A progeny-provenance trial of 6-year-old Scots pine, comprising eight native provenances each with four families in six blocks, was translocated in April 2013 to a clear-fell site in Galloway adjacent to a DNB-infected forest. Susceptibility to D. septosporum, measured as DNB severity (estimated percentage nongreen current-year needles), was assessed visually over 2 years (2013-2014 and 2014-2015; two assessments per year). There were highly significant differences in susceptibility among provenances but not among families for each annual assessment. Provenance mean susceptibility to D. septosporum was negatively and significantly associated with water-related variables at site of origin, potentially due to the evolution of low susceptibility in the host in response to high historical pathogen pressure.

Comparative transcriptomics of a complex of four European pine species.

BACKGROUND: Pinus sylvestris, P. mugo, P. uliginosa and P. uncinata are closely related but phenotypically and ecologically very distinct European pine species providing an excellent study system for analysis of the genetic basis of adaptive variation and speciation. For comparative genomic analysis of the species, transcriptome sequence was generated for 17 samples collected across the European distribution range using Illumina paired-end sequencing technology. RESULTS: De novo transcriptome assembly of a reference sample of P. sylvestris contained 40968 unigenes, of which fewer than 0.5% were identified as putative retrotransposon sequences. Based on gene annotation approaches, 19659 contigs were identified and assigned to unique genes covering a broad range of gene ontology categories. About 80% of the reads from each sample were successfully mapped to the reference transcriptome of P. sylvestris. Single nucleotide polymorphisms were identified in 22041-24096 of the unigenes providing a set of ~220-262 k SNPs identified for each species. Very similar levels of nucleotide polymorphism were observed across species (π=0.0044-0.0053) and highest pairwise nucleotide divergence (0.006) was found between P. mugo and P. sylvestris at a common set of unigenes. CONCLUSIONS: The study provides whole transcriptome sequence and a large set of SNPs to advance population and association genetic studies in pines. Our study demonstrates that transcriptome sequencing can be a very useful approach for development of novel genomic resources in species with large and complex genomes.

Rangewide ploidy variation and evolution in Acacia senegal: a north-south divide?

Knowledge of rangewide variation in DNA content and ploidy level may be valuable in understanding the evolutionary history of a species. Recent studies of Acacia senegal report diploids and occasional tetraploids in the Sudano-Sahelian region of sub-Saharan Africa, but nothing is known about the overall extent of DNA ploidy variation within the species. In this study, we determine the DNA content and ploidy level of A. senegal across its native range, and explore whether the variation is related to its evolutionary and colonization history. We used propidium iodide flow cytometry (FCM) to estimate DNA content (2C value) and infer ploidy in 157 individuals from 54 populations on various tissues, using seeds, fresh leaves, dried leaves and twigs and herbarium specimens. The mean 2C DNA (pg ± s.d.) contents detected were 1.47 ± 0.09, 2.12 ± 0.02, 2.89 ± 0.12, and a single individual with 4.51 pg, corresponding to a polyploid series of diploid, triploid, tetraploid and hexaploid individuals. Diploids were confirmed by chromosome counts (2n = 2x = 26). Most populations (90.7 %) were of single ploidy level, while mixed ploidy populations (9.3 %) comprising mostly diploids (2x+3x, 2x+4x and 2x+6x) were restricted to the Sudano-Sahelian and Indian subcontinent regions, its northern range. The species is predominantly diploid, and no mixed ploidy populations were detected in east and southern Africa, its southern range. The geographic pattern of ploidy variation in conjunction with existing phylogeographic and phylogenetic data of the species suggests that polyploids have occurred multiple times in its evolutionary and recent colonization history, including contemporary ecological timescales. The successful use of external tissues of dried twigs in FCM is new, and presents the opportunity to study numerous other dryland woody species.